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Abstract 

We propose a mechanism that can create a mass gap in the SU(2) chiral spin model at 
arbitrarily small temperatures. We give a sufficient condition for the mass gap to be non- 
zero in terms of the behaviour of an external Z(2) flux introduced by twisted boundary 
conditions. This condition in turn is transformed into an effective dual Ising model with 
an external magnetic field generated by SO (3) vortices. We show that having a nonzero 
magnetic field in the effective Ising model is sufficient for the SU(2) system to have a mass 
gap. We also show that certain vortex correlation inequalities, if satisfied, would imply a 
nonzero effective magnetic field. Finally we give some plausibility arguments and Monte 
Carlo evidence for the required correlation inequalities. 



Research supported in part by NSF Grant PHY89-15286 
2 e-mail address: kovacs@physics.ucla.edu 



1 



1 Introduction 



Two dimensional lattice spin models with a continuous non-Abelian symmetry are commonly 
believed to have no phase transition at finite temperature. According to the Mermin- Wagner 
theorem in these systems there is no ordered low temperature phase with spontaneous break- 
down of the global symmetry. However this does not in itself rule out the presence of a phase 
transition. Indeed the 2d XY model does have a finite temperature Kosterlitz-Thouless type 
phase transition characterised by power-law decay of the correlation function below the crit- 
ical temperature as opposed to an exponential fall-off in the high temperature phase. Two 
is the lower critical dimension for these models and the actual phase structure depends on 
the properties of the internal space in a very delicate way. 

In the case of spin models with Abelian symmetries there are rigorous results concerning 
their low temperature behaviour. For discrete symmetries (e.g. the Ising model) there is a 
distinct low temperature phase with long range order while for U(l) symmetry (XY model) 
the presence of the KT transition has been rigorously proved Recently the exact mass 
gap of two dimensional spin models with various different non-Abelian symmetries has been 
calculated using the Bethe ansatz 0. Scaling of the correlation function has also been 
observed in a high order strong coupling expansion of the N = oo SU(N) xSU(N) chiral 
model ||. It would however be very desirable to have a rigorous proof of a nonzero mass gap 
in these models starting from first principles without any further assumptions. This would 
probably also contribute to our intuitive understanding of the low temperature behaviour of 
non-Abelian spin models. 

The situation is possibly even worse in 4d gauge theories. While the U(l) theory has 
been rigorously proved to have a deconfining transition |4| , 4d non-Abelian gauge theories are 
believed to remain confining down to arbitrarily small couplings but a rigorous proof of this 
has not yet been found. In this case we do not even have anything analogous to the Bethe 
ansatz solution of two dimensional spin models. In view of the analogies between 2d spin 
models and 4d gauge theories it would be highly desirable to have a unified physical picture 
of their low temperature behaviour. In particular we would like to understand confinement 
in gauge theories and the nonzero mass gap in 2d spin models on a similar ground. 

For a few years there has been an ongoing programme to achieve this both in 4d gauge 
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theories || and in the analogous 2d spin models J], [7|. In the present paper we report on 
recent progress in this direction in the context of the SU(2) principal chiral model. In this 
model the degrees of freedom are SU(2) group elements attached to the sites of a square 
lattice with a ferromagnetic nearest neighbour interaction tr^l^)- We shall propose a 
mechanism that is sufficient for this system to maintain a mass gap at arbitrarily low nonzero 
temperature. 

Instead of studying the spin-spin correlation function we shall give a sufficient condition 
for the mass gap to be non-zero in terms of the behaviour of an external Z(2) flux introduced 
by twisted boundary conditions. This allows us to explicitly separate the Z(2) degrees of 
freedom belonging to the centre of SU(2) and to establish the equivalence of the SU(2) spin 
model with a dual Ising model interacting with SO (3) spins on the original lattice. We 
show that a non-zero magnetisation in the dual Ising model would imply a mass gap in the 
original SU(2) model. Due to their coupling to the Ising spins, SO (3) vortices will turn out 
to generate an external magnetic field for the Ising spins, provided certain vortex correlation 
inequalities are satisfied. As we shall see, up to this point all of our arguments are rigorous. 
Finally we shall also give some plausibility arguments and Monte Carlo evidence that the 
vortex correlation inequalities are indeed true. This gives a complete picture of the mass 
gap generation in the SU(2) spin model. 

The plan of the paper is as follows. In Section 2 we introduce the order parameter 
that we want to study, the twist. This turns out to be technically more advantageous than 
looking at the asymptotic behaviour of the spin-spin correlation function and it can be 
proved that "massive" behaviour of the twist implies the same for the correlation function. 
Also in Section 2 we shall briefly review some properties of the twist, in particular its high 
and low temperature expansion and its behaviour under duality. In Section 3 we describe 
the separation of Z(2) and SO (3) degrees of freedom and discuss the relevant SO (3) vortex 
excitations and their dynamics. By an additional Z(2) duality transformation we map the 
original SU(2) system on a dual Ising model interacting with SO (3) spins. We also rewrite 
the twist in this language. In Section 4 we develop an effective theory of the dual Ising 
system by substituting the effect of the SO (3) spins with an effective Ising coupling and a 
magnetic field. More precisely we derive sufficient conditions for the disorder correlation 
in an effective Ising model to be an upper bound of the twist in the original model. These 
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conditions turn out to be vortex correlation inequalities containing a parameter, the effective 
magnetic field of the Ising system. For zero magnetic field they are easily seen to be true 
and in Section 5 we argue that a small but nonzero magnetic field can be chosen so that the 
inequalities remain true independently of the lattice size. Also in Section 5 we present some 
Monte Carlo evidence supporting this. In Section 6 we draw our conclusions and make some 
final remarks. 

2 The twist as an order parameter 

We want to distinguish between a massive and a massless phase in the SU(2) spin model. 
By definition a massive phase is characterised by an exponential fall-off of the spin-spin 
correlation function (|tr UqU x ) at large distances x, where () means expectation in the 
infinite volume limit. In our framework it turns out to be technically more adventageous 
to consider another operator, a twist winding around the lattice and its behaviour as a 
function of the finite lattice size.0 Since the twist has not been very commonly used as an 
order parameter in spin models (but see e.g. 0), in this section we shall collect some useful 
results about it to make the paper self-contained. 

2.1 Notation 

Although we use notations that have become more or less standard in (lattice) field theory, 
we find it useful to include here a brief section on the notation especially for later reference. 

We shall work on a finite two dimensional periodic square lattice A. Sites, links and 
plaquettes of the lattice will be denoted by s, I and p. Because of its particular simplicity 
we shall make use of the language of group valued forms when writing down the fields and 
their interactions. By a G-valued n-form we mean an assignment of G elements to oriented 
elementary n-dimensional simplices of the lattice. Thus e.g. a Z(2) valued 2-form u assigns 
a Z(2) element to each plaquette and a Z(2) element on a particular plaquette p will be 
denoted by uj v . 

3 The analogous quantity in gauge theories is the sourceless 't Hooft loop or magnetic flux free energy 
that has been extensively used to describe the phases of lattice gauge theories § . 
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There are two useful operations that can be defined on forms, the exterior derivative "d" 
and its dual "d*". They map an n-form on an n — 1-form (d) and n + 1-form (d*) respectively. 
By definition 

{du) k = \{ u m , (1) 

where k denotes a link, m a site if 00 is a 0-form; k is a plaquette and m is a link if u is a 
1-form and d is the boundary operator. Similarly 

(dV) m = n ^, (2) 

k:m£dk 

where m is a site, A; is a link if uj is a 1-form and m is a link, A; is a plaquette if a; is a 2-form. 
Thus e.g. the curvature of a Z(2) valued gauge field, a 1-form A, can simply be written as 
dA. For better readability we shall always omit the parentheses so e.g. du m means (dcu) m . 

In most of the cases we shall only use Z(2) valued forms the only exception being the 
dynamical variable in the SU(2) chiral model which is an SU(2) valued 0-form U s (spins live 
on lattice sites). The interaction between the two spins residing at the ends, si and s2 of a 
given link I, can be written as dUi = tr(Ul 1 U s2 ). This is the only case when the orientation 
of elementary simplices will matter. 

When calculating thermal averages involving (discrete or continuous) group valued fields 
integration is always understood with respect to the normalised Haar measure. 

2.2 Twist and duality in the Ising model 

At first we look at the twist in the context of the Ising model since this will be needed later 
when we discuss the effective Ising model of the SU(2) spin system. In the Ising model the 
twist is not particularly useful for distinguishing between the phases of the theory since the 
magnetisation provides a much simpler order parameter. Nevertheless as we shall see, the 
twist can also be used for this purpose. 

The partition function of the Ising model is 

Z = Y[Jd<T s exp(p'E i d<T?J (3) 

with spins a s — ±1 living on a finite periodic square lattice. Here / da s means integration 
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with the normalised invariant measure?] on Z(2), sites and links of the lattice are labelled by 
s and I and finally dai is the product of the two spins on the link I. 

By definition the twist r% — ±1 along the direction i is the operator that changes the 
couplings from (3 to Ti/3 on a stack of links winding once around the lattice along the given 
direction (Fig [l]). Physically r» = — 1 means that a topologically nontrivial "do main- wall" 
was created along the affected links. Although this domain-wall is closed, it is not a boundary 
of any region therefore it is impossible to transform it away by a suitable redefinition of some 
of the spin variables. In contrast, if the domain- wall were a boundary of a region then by 
a change of variables a — > Ti<j in the region bounded by the twist, it would be possible to 
cancel it and the partition sum would not depend on Tj. 

By a similar argument it can be easily seen that the twisted stack of links can be "con- 
tinuously" deformed by a change of variables and the twisted partition sum does not depend 
on the actual location of the twisted links as long as they form a closed loop winding once 
around the lattice. 

The two phases of the system can be characterised by its response to such a twist. 
Intuitively one expects that in the high temperature disordered phase, having a twist does 
not cost too much energy on a sufficiently large lattice. On the other hand in the low 
temperature broken phase where most of the spins tend to be in one direction, the free 
energy of the twist is expected to grow with the lattice size. 

To make these ideas more precise we can consider the following observables 

G{L) - 2 Z + (L) + Z_(L) and G{L) ~ l(Z.(L) + Z + (L)Y (4) 

Z ± {L) = j ± dr 2 Z(±l,r 2 ), (5) 

where Z(t 1 ,t 2 ) is the partition function with twists t\,t 2 = ±1 in the two independent 
directions and L is the linear size of the lattice. (The log of) G*(L) measures how the free 
energy of a twist depends on the lattice size and G(L) is the Z(2) "Fourier-transform" of 
G*(L). The form of the normalisation Z_(L) + Z + (L) in the denominator was chosen for 
further convenience, we could as well have used just Z + (L). 

4 Later we shall also consider continuous groups and to make the notations uniform we use the same 
language for discrete groups too. 
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Figure 1: A twist winding around the lattice. The links on which the sign of the coupling has been 
reversed are shown. 



At high temperature (small 0) G(L) decays exponentially with increasing lattice size. 
This can be most easily seen in the high temperature expansion. Recall that the high 
temperature expansion of the Ising model is the sum of closed graphs built of links of the 
lattice. The weight of each graph is proportional to (tanh/3) # of links . In the presence of a 
twist there is an additional minus sign for each twisted link contained in the graph. Since 
the twist winds around the lattice, all small loops will intersect it an even number of times, 
so these will not contribute to G(L). The lowest order contribution comes from a loop going 
all the way around the lattice in the direction perpendicular to the twist. The shortest such 
loop contains L links and is thus proportional to (tanh / 3) i . If G(L) goes exponentially to 
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zero it means that G*(L) will go to 1 in the L — > oo limit. This is indeed what we expected, 
in the disordered phase the free energy of the twist remains bounded as the lattice size goes 
to infinity. 

In the low temperature phase with the Z(2) symmetry broken the behaviour of G(L) and 
G*(L) is different. Here the free energy cost of having a twist grows linearly with the lattice 
size and G*(L) goes exponentially to zero while G(L) goes to a nonzero constant for L large. 
This is certainly true for the ground state and it can also be proved in higher orders of the 
low temperature expansion. 

Notice that the behaviour of G(L) in the symmetric phase is similar to that of G*(L) in 
the broken phase and vice versa. This is not an accident since G and G* are dual to one 
another in the following sense. The high temperature expansion of the Ising model with 
inverse temperature (3 can be identified with the low temperature expansion of a dual Ising 
model with spins living on the sites of the dual lattice i.e. plaquettes of the original lattice 
at inverse temperature (3*. This gives a mapping between the high and the low temperature 
regimes. It can be proved that the high temperature expansion of G(L) on the original 
lattice is identical to the low temperature expansion of G*(L) defined in the dual system 
which means that G{L)p = G*(L)p*. 

Finally we note that in the presence of an external magnetic field the qualitative behaviour 
of G and G* is the same as in the low temperature phase. For sufficiently high temperature 
the high temperature expansion is convergent and the leading order exponential fall-off of 
G*(L) can be rigorously verified. 



2.3 Twist in the SU(2) spin model 

It is straightforward to generalise the Z(2) twist for any other spin model with a global 
symmetry containing Z(2). In principle we could also use any element of the global symmetry 
group to define the twist, in spin models there is no constraint (not like in gauge theories 



T0[| ) restricting the twist to lie in the centre of the symmetry group. In this paper however 



we shall consider only Z(2) twists. 
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The SU(2) xSU(2) chiral spin model is denned by the partition function 

Z= n [dU s exp [P^duA , (6) 

seA J \ leA J 

where the degrees of freedom U are SU(2) group elements on the sites of the finite periodic 
lattice A and dUi = tr(C/J[/ s /) with [ss'} being the boundary of the link I. The twists and 
the quantities G(L) and G*(L) can be defined completely analogously to the Ising model. 

The high temperature behaviour of these observables is qualitatively the same as in the 
Ising model. At low temperature however there is an important difference. In the ground 
state of the twisted Ising system all the spins are aligned and the energy of the twist is 
concentrated along a stack of links winding around the lattice parallel to the twist. On the 
other hand if the spins can take their values in a continuous manifold, the energy of the twist 
can be spread along the direction perpendicular to it. 

In this case the energetically most favourable configuration has dUi < along the twist 
with the spins roughly antiparallel. The direction of the spins along the spreading direction 
changes gradually making half of a complete turn around the lattice (Fig|2|). In the ground 
state neighbouring spins in the spreading direction make an angle proportional to 1/L cost- 
ing an energy ~ 1/L 2 and the total energy of this configuration goes as ~ L 2 x 1/L 2 on 
large lattices. We can see that the energy of a twist does not diverge with the lattice size, 
at least in the classical (zero temperature) approximation. Intuitively we expect that at 
nonzero temperature the twist costs even less since there is more disorder. This mechanism 
is analogous to flux spreading in gauge theories |Tl[| . 

We saw that both in the Ising model and in the SU(2) spin model the behaviour of G(L) 
as a function of the lattice size is similar to that of the infinite volume limit of the correlation 
function. This motivates our choice of G(L) as an indicator of the massive phase. This can 
be rigorously justified by proving that in the SU(2) model exponential fall-off of G(L) with 
the lattice size implies the same asymptotic behaviour of the infinite volume limit correlation 
function. (The proof of the analogous statement in gauge theory || can be easily adapted 
to the present casa. See fll"2"fl .) 
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twisted links 



Figure 2: The "semiclassical" ground state configuration of a twisted spin model with a continuous 
spin manifold. The energy of the twist can be spread in the direction perpendicular to the twist. 

3 Separation of the Z(2) and SO (3) variables 

In the previous section we saw that it is enough to prove the exponential decay of the twist 
G(L) as a function of the lattice size in order to establish a mass gap in the SU(2) spin 
model. Since G(L) is denned in terms of a Z(2) twist it will be helpful to disentangle the 
Z(2) degrees of freedom belonging to the centre of SU(2) from the rest of the system. In this 
Section we shall perform this task by giving a set of exact transformations that eventually 
map the SU(2) model on a dual Ising model interacting with SO(3) spins on the original 
lattice. We mainly follow but in addition we also give a more detailed description of 
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the underlying dynamics of the string and vortex excitations of the SO (3) system and their 
coupling to the Ising spins. 

3.1 Gauging the Z(2) symmetry 

As a first step we introduce new Z(2) valued link variables Ai to gauge the Z(2) part of the 
global symmetry. To maintain the equivalence of the gauged model with the original one the 
gauge field Ai has to be constrained to be pure gauge by inserting a delta function 5(dA p ) 
for each plaquatte. By definition dA p is the curvature of Ai, the product of link variables 
around the plaquette p and 5(1) = 2, 8(—l) = 0. This is however not enough to make 
Ai pure gauge. In addition we also have to constarain the two independent topologically 
nontrivial holonomies Ad = Hied to be 1. Here C\ and Ci are two loops built of links 
winding around the lattice in the two independent directions. 

Putting this all together the (nontwisted) partition function reads as 




In the second line above we used the notation 

rji = sign dL^ = sign trU^U h (8) 

and the change of variables 07 = Affi to absorb the sign of dUi into a redefined gauge field 
on each link. The variables T)d and ctq denote the product of 77's and cr's around Cj. 

The twisted partition sums Z(t\, t%) can also be rewritten in the same fashion by absorb- 
ing the extra signs coming from the twist also in the cr's choosing <Ji = A^Ti on the twisted 
links. The only difference compared to (|7|) will be that the holonomies around C\ and C2 
will pick up the twist T2 and T\ along the direction perpendicular to them. Acoordingly the 
arguments of the corresponding delta functions will be multiplied by T\ for C2 and r 2 for G\. 
Contractible loops, in particular the boundaries of plaquettes always intersect the twists an 
even number of times therefore the plaquette delta functions do not have to be modified. 
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It is now easy to compute the order parameter G(L) directly from its definition 
Jdndr 2 Z(t 1 ,t 2 ) 

= \ n fdUs n fai n^dw 1 ) ^- c \ e XP f^^w j > w 

z seA J zeA peA \ «gA / 

where Z is exactly the same as the numerator of the r.h.s. except the factor crc 2 r ]cl * s missing. 
This factor in the numerator comes from the property of the delta function on Z(2) that 



Jdr t5(tct) = a. (10) 



The advantage of the normalisation that we chose for G(L) is that in this way by integrating 
out the twist we could get rid of the delta functions constraining the holonomies around C\ 
and C 2 in the "partition function" Z . This effectively means that we integrated out the Z(2) 
part of the boundary conditions. 

The remarkable property of this form of the partition function and G(L) is that the in- 
tegrand depends on the SU(2) variables only through the SU(2)/Z(2) cosets. In other words 
it has a U s — > —U s local Z(2) gauge symmetry and the spins can be regarded SU(2)/Z(2) = 
SO (3) variables rather than SU(2) ones. Notice that the r/j's themselves are not gauge invari- 
ant but products of them around any closed loop, in particular around plaquette boundaries 
and Cj are gauge invariant. The price that we pay for this extra gauge symmetry is the 
appearance of the new Z(2) link variables <j\. 

3.2 Z(2) strings and vortices 

At this point it is in order to discuss the physical meaning of the degrees of freedom and 
excitations of the system that we obtained by rewriting the SU(2) model in terms of Z(2) 
and SO (3) variables. 

In (|9D there are two different types of couplings between the SO (3) spins U s and the Z(2) 
link variables <j\. At first <j\ contributes an extra sign to the coupling between neighbouring 
SO (3) spins residing on the two ends of the link I. If o\ = 1 then these two spins are more 
likely to be parallel whereas for o\ = —1 they tend to be "perpendicular". In addition the 
plaquette delta functions constrain the "curvature" da p = Yiiedp °i to be equal to drj p . 
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The relevant Z(2) excitations of the model are (stacks of) links with o\ = —1, which 
we call cr-strings. These strings are either closed or they terminate on plaquettes having 
dcr p = — 1 i.e. an odd number of negative cr's on their boundary. Due to the plaquette delta 
functions the end-plaquettes of cr-strings have to coincide with the df] p = — 1 plaquettes. For 
reasons that will become clear later we shall call these vortices. 

In analogy with cr-strings we can also define r/-strings as the location of rji = —1 links. 
By construction every vortex is an endpoint of both a cr and an r/-string (Fig |3[). Notice 
however that the location of the //-strings does not have a physical meaning in terms of the 
SO (3) variables, they can be deformed by U s — > —U s gauge transformations. The location 
of vortices on the other hand is gauge invariant since Z(2) gauge transformations always 
change the sign of an even number of factors in dr] p . The situation is analogous to gauge 
theories where one has Dirac strings attached to monopoles and the strings can be deformed 
by gauge transformations but their endpoints (boundaries in 4d), the monopoles (monopole 
loops) are gauge invariant objects. 

It should be noted that there is a substantial difference between the energetics of the 
a and rj strings. To see this let us look at a closed contour of links enclosing exactly one 
vortex. Both the cr and the r] string attached to the vortex have to pierce the contour 
somewhere. The link I where the a string intersects the contour is unambiguously given by 
the [a] configuration and it carries an energy ~ 2 dlli. On the other hand as we have already 
seen, the location I' where the r/ string crosses the contour is not Z(2) gauge invariant. It is 
thus not surprising that I' has no distinguished role among the links of the contour; it does 
not carry any extra energy. Indeed in a typical low temperature SO (3) configuration the 
spins vary slowly around the contour and they make half of a compelete turn as the contour 
encircles the vortex. This is more favourable both in terms of energy and entropy than 
having abrupt changes. Now we can see why the dr/ p = — 1 plaquettes were called vortices. 
dr] p is well defined in terms of the SO (3) variables and a "smooth" vortex configuration is a 
special case of defects appearing in 2d spin models due to the non simply-connectedness of 
the internal space [TT3 . 



At this point we want to emphasize that the internal space of the original SU(2) spin 
model i.e. the SU(2) group manifold is simply connected, therefore no vortices are present 
in the original model. After separating the Z(2) degrees of freedom however, the remaining 
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® vortex o\ = —1 r\\ = — 1 

Figure 3: A typical configuration with one vortex pair and the r\ and a string connecting them. 

configuration space (at each site) is essentially SO (3). Due to the doubly connected nature 
of SO (3) this representation admits vortices characterised by a Z(2) charge. There is no 
contradiction between the two representations of the model because the SO (3) vortices always 
come along with da p = — 1 defects in the Z(2) part of the system (this is ensured by the 
constraints in eq. (|7|)). Were we to rewrite the model again in the original SU(2) language 
the two types of defects, the vortices and the Z(2) defects in a, would always combine to 
give defects in the SU(2) language that are topologically trival. 

Let us now consider two vortices connected by an r\ and a a string. It is clear that at low 
temperature the interaction between the two vortices is dominated by the a string since its 
energy is proportional to its length. On the other hand from a simple semiclassical estimate 
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(i.e. using a slowly varying "smooth" SO (3) configuration) the energy of the 77 string can be 
easily seen to go as the log of its length. At low temperatures the vortices are in closely bound 
pairs due to the linearly confining potential of the a string. Vortex pairs with short a strings 
between them are local excitations, they are not expected to have a big influence on the large 
scale behaviour of the system. The 77 strings however can still fluctuate considerably because 
of their smaller energy cost and as we shall see long 77 strings are capable of disordering the 
system even on large distance scales. 

At first sight one would expect that at sufficiently low temperature due to their loga- 
rithmic behaviour even the long 77 strings will freeze out of the system. This is however not 
quite right. To see this we have to notice that in the above semiclassical energy estimate 
we assumed that the lattice is very large compared to the length of the 7/ string. Since we 
shall not be working directly in the infinite volume limit but rather consider the lattice size 
dependence of G(L), this assumption may not be right. Indeed, let us look at two nearby 
vortices with a short a string in between but the 77 string connecting them going all the way 



around the lattice. By a semiclassical argument similar to the one in Section |2.3| we can see 
that the energy of this configuration goes to a constant when L — > 00. The two vortices of 
course can be removed by closing the r\ string around the lattice. 

These types of configurations can potentially have a long range disordering effect. In 
fact it is exactly these configurations that have to be present with a sufficient weight to 
ensure that the relative difference between the twisted partition sum and the untwisted one 
decreases sufficiently rapidly with the lattice becoming larger. As we shall see later, this is 
what we need for G(L) to decay exponentially. 



3.3 Duality transformation 

In the previous section we have succeded in separating Z(2) and SO (3) variables in the 
SU(2) spin model and showed that vortex and string excitations of the SO (3) part can have 
a long range disordering effect. Unfortunately the Z(2) part of the system does not look very 
familiar and in order to proceed we want to cast it into a more managable form. This can 
be done by a duality transformation on the a link variables which amounts to trading them 
for Z(2) plaquette variables that we call uo p . The interaction between uu's will turn out to 
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be nearest neighbour ferromagnetic although with fluctuating couplings that depend on the 
SO (3) spin configuration. The a/s in this way give an Ising model on the dual lattice the 
sites of which are plaquettes of the original lattice. 

Technically the Z(2) duality transformation is done by expanding each factor in the 
integrand of G(L) (equation @) depending on the cr's in terms of Z(2) characters as 



f(a) = /(l) + af(-l) = J da f(a) Xa (a). (11) 
The expansion for the different factors in (||) reads as 

^(dopdr^ 1 ) = Jdujp x^idcrpdr)- 1 ) 



e 



&c 2 = X-i{vc 2 ) = II X-lfa 

°i\dU,\ _ 



da lXai (vi)e cm > a '\ (12) 



where e c ^ dUl ' ai ^ is the Z(2) "Fourier transform" of e CTi ' dl/i ' and for later convenience it can be 
split into a part depending only on the SO (3) variables and another one depending both on 
the SO (3) and Z(2) degrees of freedom; 

e C{AU u ai) _ 1 ^|dDi| + ai e -(3\dUi\^ = e Af(dL/,) e «,Jr(dL/,). ( 13 ) 

Here the functions K and M are given by 



M(dU[) = -lnednh(2^|dZ7i|) and K{dU{) = - lncoth(/3|dt/i|). (14) 

Substituting these into the expression of G(L) the a variables can be explicitly integrated 
out using the orthogonality of the characters. This will give rise to constraints between the 
remaining Z(2) degrees of freedom, namely on each link 

ai = d*ui = Yl u p (15) 

except on the links belonging to C<i where ai = — d*Ui. These constraints make the ai 
integrals trivial and yield 



G(L) = i [du[U] n JduJpX.Mn, 



x Vc2 exp [ ]T K(dU t ) d*ujt - K(dU t ) d*u t ] , (16) 
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where 

du[U] = H^.^S^ (17) 

seA 

is the SO (3) (Z(2) gauge invariant) part of the measure. 

Exactly the same transformation can be carried out on the partition function Z, the only 
difference in the final result compared to (|16D will be the absence of the factor t]c 2 and the 
minus sign in the couplings between cj's along Ci- 

In this representation the system consists of Z(2) spins u p attached to plaquettes and the 
SO (3) part of the system is unchanged. The cj's on nearest neighbour plaquettes sharing the 
link I interact via the coupling K(dUi). The Z(2) part of the system is essentially an Ising 
model on the dual lattice although the spin-spin couplings can fluctuate since they depend 
on the SO (3) configuration. 

For large (3 the effective SO(3) action M(dUi) is peaked at the maximum of \dlli\ which 
means that /3|dL^| is also large for "most" of the relevant configurations. The asymptotic 
behaviour of M{dU{) in this limit is ~ /3|dC//| — In 2/2. On the other hand the duality 
transformation changes low temperature to high temperature for the Z(2) part of the system. 
Indeed the (3 — > oo limit of the effective Ising coupling is K(dUi) ~ e~ 2/3 ' d[/i '. 



4 Construction of the effective Ising model 

We have seen that after separating the Z(2) and SO(3) variables in the SU(2) model, G(L) 
can be rewritten in terms of a dual Ising model coupled to SO (3) spins. As can be seen from 
equation ([Tj]) the Ising part of G(L) is essentially the ratio between a twisted (along C2) and 
the untwisted partition sum, i.e. the free energy of a twist along C 2 (G*(L) in the language 
of Section [|). Notice that because of the duality transformation, G(L) in the original model 
is becomes G*(L) in the dual Ising model. 

We have already discussed that the Ising couplings between nearest neighbour plaquettes 
depend on the SO (3) spin configuration. Besides there are two more SO (3) dependent pieces 
in the expression of G(L); T]c 2 and the product of group characters Xu> p (dri P ). While i]c 2 
depends solely on the SO (3) variables, the group characters couple vortices to the Ising 
spins on each plaquette. 
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For a fixed SO(3) configuration with vortices at (pi,P2---P2n) an d S r\ strings^ crossing 
G% these give an overall factor of (— l) s Yli=i u Pi- Were it not for this additional factor 
depending on the SO (3) configuration, the Z(2) part of the system would be a ferromagnetic 
Ising model at high temperature with unbroken symmetry and G(L), the ratio of the twisted 
and untwisted partition sum would go to a nonzero constant on large lattices. It follows that 
if we were to eliminate all the vortices from the measure and also constrain rjc 2 to be +1 
the exponentially falling asymptotic behaviour of G(L) would be lost. The configurations 
responsible for the correct asymptotic behaviour of G(L) are exactly those that would be 
eliminated by the above constraints. 

4.1 Strings and vortices in the dual representation 

Let us now have a closer look at the relevant configurations in the different vortex sectors. 

In the absence of vortices the only difference between the numerator and the denominator of 

fll6|) is the factor rjc 2 - This is true up to graphs of size ~ L in the high temperature expansion 

of the us Ising model but these are exponentially suppressed (see Section |27^ ). All the short 

contractibe rj strings have to cross an even number of times and they do not contribute 

to r]c 2 - f]c 2 = ~ 1 thus means that there is an odd number of r\ strings winding around the 

lattice perpendicularly to C%- In the zero- vortex sector it is exactly the long topologically 

nontrivial 77 strings that give a negative contribution to the numerator of G(L) and a positive 

one to the denominator. If there were no vortices at all, the asymptotic behaviour of G(L) 

would solely depend on the relative weigth of these configurations. 

Let us now consider the sector with two nearby vortices at p\ and pi- In this sector the 

high-temperature expansion of the u spin system contains small graphs connecting p\ and 

Pi- If a graph like this crosses C2 then in the numerator it acquires a minus sign due to the 

twist. Besides we still have the factor r] C2 . The combined effect of these two signs will be 

different for the numerator and the denominator of ( |l6l) only if the uj graph together with 

the 77 string connecting the two vortices close into a noncontractible loop (Fig This is 

absolutely essential, otherwise the sign coming from r]c 2 would cancel the one coming from 

D Although S is not well denned in terms of the SO(3) variables, its parity (— l) s is invariant under coset 
reparametrisations 
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(g) ^2 

<8> 



(g) vortex rji = — 1 

cj- string 

Figure 4: A vortex pair with its connecting r] and u string forming a noncontractible loop winding 
around the lattice. This configuration gives different contributions to the numerator and denomi- 
nator of G(L). 

the twist. In other words the curve C2 could be deformed by a Z(2) gauge transformation 
to decouple from both the u and the r\ string (as in Fig [5]). This obviously cannot be done 
with the configuration in Fig [|. 

The configuration in Fig [| is not very much different from the one discussed in the 0- 
vortex sector. Again it contains a long rj string winding almost around the lattice connecting 
the two vortices with a small gap between them. 

We can already see the similarity between this expansion and the high temperature 
expansion of the Ising model with a nonzero magnetic field discussed in Section [2.2| . The 
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<8> 




c 2 


<8> 







(g) vortex rji = — 1 



cj-string 

Figure 5: A vortex pair with the connecting r\ and u string forming a contractible loop. This 
graph gives the same contribution to the numerator and the denominator of G(L). 

only difference here is that the graphs of order In are proportional to In- vortex expectations 
instead of the 2n-th power of the magnetic field. It should be stressed that the vortex 
expectations have to be calculated with the SO (3) measure (|T7| ) which is "blind" to the 
centre of SU(2). 

The similarity between the two expansions hints that it might be possible to approximate 
the dual representation of the SU(2) model with an effective Ising model. In this way we could 
trade the SO (3) variables for a suitably chosen effective coupling K e ^{(3) and an effective 
magnetic field h(/3). To put it more precisely, we want to bound expression ([16]) of G(L) by 
the ratio of the twisted and untwisted partition function of an Ising model with a suitably 
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chosen coupling K e s(/3) and a magnetic field h{(3). For large f3 K e s and h will be exponentially 
small and the high temperature expansion for the effective Ising model will be convergent. 
As we saw in Section |2]^ in the presence of a nonzero magnetic field the ratio of the twisted 
partition function to the untwisted one falls off exponentially with the lattice size. It means 
that if the Ising bound for G(L) can be established with a nonzero magnetic field it would 
imply the exponential decay of G(L). 

To summarize, in order to explicitly construct the effective Ising model we have to deal 
with three questions. 

• The fluctuating couplings K{dUi) have to be approximated by some K e{i . 

• The effect of the SO (3) vortices on the Ising spins has to be substituted with an external 
magnetic field. 

• Somehow we have to get rid of the factor 7]c 2 that depends only on the SO (3) variables. 

The last point turns out to be the simplest one. Since long r\ strings disorder the system, 
constraining t]c 2 to be +1 makes the system more ordered and increases G(L). Indeed by 
using reflection positivity it can be proved that this constraint increases G(L) |T2[ . From 
this point on we shall consider the quantity G(L, C^ 1 ") which is defined as G(L) calculated 
with this additional constraint. The SO(3) measure dv\U\ supplemented with this constraint 
will be called dv + [U}. 

4.2 The effective Ising coupling 

We now want to compare G(L, C^) with the quantity 




i.e. the ratio of the twisted and untwisted partition function of an Ising model. 

In the original SU(2) model we started with the particular form of the interaction dUi = 
trlll^i- We could have as well used any other action in the same universality class. In 
particular the choice 

leA zeA 
21 



with fi > would have considerably simplified the discussion. After the separation of the 
Z(2) variables and the duality transformation on this action we get 

M(dU t ) = p\dUt\ and K(dU t ) = - lncoth 2///3, (20) 

i.e. the Ising coupling is independent of the SO (3) variables and we do not need to find an 
approximation for K e g. Because of its simplicity we shall only consider this form of the 
action and the quantity G{C^ ft, L) corresponding to G(L, C%) with this particular choice 
of the SU(2) action. 



4.3 Ginibre decomposition 

We now want to find a nonzero magnetic field h, independent of the lattice size such that 
on any lattice 

G cS (h, K, L) - G(C+, ft, L) > 0. (21) 

This would imply the exponential fall-off of G(C , 2 I ", /x, L) and by the inequalities proved so 
far the same behaviour for the correlation function of the SU(2) spin model. 

By an argument similar to Ginibre's proof of the Griffiths' inequalities JM the difference 



Z e gZ(G e f[ — G) can be decomposed into a sum of terms independent of the u variables, 
containing only vortex correlations calculated with the SO (3) measure du + [U], including the 
constraint on the rj string. Details of this decomposition can be found in the Appendix. The 
form of these terms is (up to non-negative constant factors) 

D(h) = (i[(e; ± he;) n_(e; ± he;)) L , (22) 

where ()l means expectation with respect to the SO(3) measure dv+[U], A is some subset 
of A containing an even number of plaquettes and h = tanh h. 

It is easily seen that for h = all these diagrams are strictly positive implying that 
G e ff(h = 0, K, L) > G{C/2 i fi, L). On any finite lattice G e g is an analytic function of h and 
this inequality is in fact true not only for h = but for any h in some finite interval around 
0. However we still need to establish that this interval does not shrink to zero as the lattice 
size goes to infinity. 
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5 Vortex correlations 



In this section we want to study the dependence of the vortex expectations ( p2|) on the 
effective magnetic field and make it plausible that there is a finite interval around h = for 
which all these expectations are non-negative independently of the lattice size. We would 
like to emphasize that up to this point all our arguments were rigorous. 



5.1 Factorisation inequality 



Vortex expectations of the type ( p2|) come in huge varieties. We can introduce some order 
into this abundance by explicitly constraining out some of the Z(2) excitations while still 
keeping a finite density of them. This will certainly make the system more ordered and the 
presence of a mass gap in this system would imply a mass gap in the original model. With 
these further constraints we can essentially pair up vortices so that the remaining nonzero 
expectations of the type (^2|) will have the form 

D(h) = (\[{e Pi -he + Pi ) n^ P + ), (23) 

Pi£A p^A 

where A is a sublattice of A, Pj = {pnPi2} denotes pairs of plaquettes and 9$. = 9p a 0p i2 - 
Now for h = it can be proved that 

n 

<n ^lk> > n<^ ik + >- (24) 

P,GA pf^A 1=1 viPi 

This correlation inequality can be checked by using reflection positivity for vortex configu- 
rations that are symmetrically placed about a line bisecting the lattice. For more general 
configurations the proof involves the application of the FKG inequalities . 



Now it is quite plausible that if the above factorisation inequality holds for h = then 
it will be true for some finite interval around h = which is independent of the lattice size. 
Let us for the moment assume this. It follows then that we have a lower bound on G e g — G 
in terms of a sum that contains products of two-vortex expectations of the form 

n 

i=i p^p, 
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It is now enough to establish that all these factors are non-negative for any h in some finite 



interval around zero independent of the lattice size. This is equivalent to the statement that 
the free energy of a vortex pair (as compared to that of the no- vortex state), 



is bounded for arbitrarily large lattices. Normally one would think about a pair of vortices 
as a local excitation and its free energy is not expected to diverge with the lattice size. This 
is however not trivially true in our case due to the constraint on the r\ sting contained in 
the measure <iz/ + [[/]. If the two nearby vortices are cut off from each other by C2, their 
connecting rj string has to wind around the lattice to avoid crossing C<i- This type of 
excitation is definitely not a local object. In fact the "most non-local" excitation in the 
two- vortex sector is the one that has two adjecent vortices separated by C2 (see Figure f|). 
All other types of two-vortex configurations contain shorter r] strings and therefore have a 
smaller free energy. 

In two dimensions, as we have already seen, the free energy of such a long 77 string stays 
finite as L — > 00 in the semiclassical approximation. This is due to flux spreading that allows 
the cost of the string creation to spread laterally in the direction perpendicular to the string. 
Intuitively one expects the spreading of the flux to be even faster than that given by the 
semiclassical approximation but unfortunately we could not prove this analytically. Instead 
we measured the lattice size dependence of the two- vortex free energy using Monte Carlo. 

5.2 Monte Carlo results 

In this subsection we want to give some evidence that the free energy of two adjecent vortices 
separated by the curve C2 does not diverge when the lattice size goes to infinity. Moreover 
we claim that this is very likely to be true for any non-zero temperature. 

It is well-known that both in non-Abelian 2D spin models and the analogous 4D gauge 
theories the strong and weak coupling regions are separated by a crossover region, where 
the specific heat has a finite peak |TB|. We used the specific heat peak to ensure that the 
coupling that we use in the simulations is already "weak". The Monte Carlo simulation 
was performed at j3 = 2.0 which is already in the weak coupling region as can be seen from 




(26) 
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inverse temperature (beta) 

Figure 6: The specific heat and the vortex density versus inverse temperature (/?) on an 8x8 lattice 
calculated with the SO(3) action. 



Figure || We also plotted in the same Figure the Z(2) vortex density versus the inverse 
temperature. Another characteristic feature of the weak coupling region is that the vortex 
density (monopole density in gauge theories) decreases exponentially as a function of (3 f!6 
This can also be clearly seen in the Figure. 



We do not expect any dramatic phenomenon to happen when going from the weak cou- 
pling side of the specific heat peak to even weaker and weaker couplings, therefore our results 
should qualitatively apply for any "weak" but of course non-zero coupling. 

We measured the probability of the configuration with two adjecent vortices at fixed 
positions separated by C2, normalised by the probability of having no vortices at all. The 
calculation was done by generating a series of configurations using a local heat bath algorithm 
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with the S0(3) invariant measure g?i/[?7] but without the constraint on the r\ string. In a long 
Monte Carlo run the number of configurations of the type in Figure |] was counted and 
divided by the number of configurations with no vortices at all and an even number of 77 
strings crossing Ci- 

Since (at fixed 0) on larger lattices there are typically more vortices, the probability of the 
configurations with exactly zero and one vortex pair decreased very rapidly with increasing 
lattice size. This meant that the quantity we wanted to measure was given as a ratio between 
two numbers both becoming very small on larger lattices, however their ratio was expected 
to be stable. This made the signal less accurate on larger lattices. An improvement by a 
factor L 2 could be achieved by counting all the configurations that were translations of the 
one in Figure [| and dividing the result by L 2 . Because of the periodic boundary conditions 
the lattice had an exact translation invariance and this procedure did not change the results. 

The simulations were performed on square lattices of 5 < L < 13. The rapidly deterio- 
rating quality of the signal made it impossible to go beyond L = 13. Even at this point we 
typically needed several hundred thousands of independent configurations to get a signal at 
all. Our results are summarised in Figure 0. We can see that the probability of a vortex 
pair with a long 77 string decreases on small lattices until it stabilises at a moderate lattice 
size (L ~ 8 — 9) and stays constant thereafter. Recall that for our purposes it is enough if 
this quantity remains non-zero in the L — > 00 limit. 

To study the spreading of the flux of the long eta string, we measured the same quantity 
with the lattice size fixed in the direction parallel to the string and varied in the perpendicular 
direction. The results are plotted in Figure ||. The pronounced effect of flux spreding is 
obvious; as the lattice becomes wider, there is more space available for the flux of the r) 
string to spread and its probability increases dramatically. 

To summarise, our Monte Carlo results are qualitatively consistent with the semiclassical 
picture of flux spreading and it is indeed very hard to imagine that anything could happen 
either at higher (3 or larger lattice sizes that could make the two-vortex probability vanish 
in the L — > 00 limit. Of course our results do not imply any particular analytic form as to 
how the flux actually spreads and it seems very hard to distinguish between a semiclassical 
"massless" spreading and an exponential spreading. This is however not necessary for our 
purposes; the only property that we need is that the probability of the long string with two 
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lattice size 

Figure 7: The probability of a neighbouring vortex pair with its r] string winding around the lattice 
(see Fig |j) as a function of the lattice size at f3 = 2.0. The measurement was done with the SO(3) 
action. 



vortices remain non-zero in the L — * oo limit. 



6 Conclusions 

We proposed a mechanism that is sufficient to create a nonzero mass gap in the SU(2) chiral 
spin model at arbitrarily small temperatures. Instead of considering the correlation function 
we were looking at an operator (G(L), the twist) that measured how effective fluctuations 
were in destroying the correlations between the "relative signs" of spins at different locations. 

In this way we could separate the spin system into two interacting parts, a Z(2) and an 
SO (3). We derived the conditions that the SO (3) part had to satisfy so that the mass gap 
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Figure 8: The same as in Figure [/] except that instead of using a square lattice the length in the 
direction parallel to the r/ string was fixed to 9 and only the transverse (perpendicular to the string) 
size was varied. 



could be rigorously established in terms of an effective Z(2) system, which was an Ising model 
on the dual lattice. By using a plausible (but so far not proved) vortex correlation inequality 
we could reduce these necessary conditions to just one condition, namely the boundedness 
of the free energy of a vortex pair connected by an rj string winding around the lattice for 
any lattice size. 

This condition on the truly non-Abelian part of the system was substantially weaker 
than the existence of a mass gap itself. The presence of a mass gap would require that an 
external Z(2) flux, introduced by twisted boundary conditions, spread exponentially fast on 
large lattices, i.e. the expectation of a twist would have to go to 1 exponentially. On the 
other hand in our scheme we only needed that the expectation of the twist be non-zero for 
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asymptotically large lattices. It is nevertheless surprising that we failed to verify even this 
substantially weaker condition rigorously and eventually we had to resort to Monte Carlo 
to check it. As expected, our Monte Carlo data show that even at low temperature the 
expectation of the twist goes to a nonzero constant as L — > oo. It would be however very 
desirable to have an analytic proof of this and the vortex correlation inequalities too. 

This framework can also be applied to 3d and 4d SU(2) gauge theories with essentially no 
modifications. The only minor difference here is that the corresponding objects on the lattice 
live on higher dimensional simplices. For example instead of Z(2) vortices on plaquettes we 
have Z(2) monopoles (monopole loops in 4d) on cubes and the twist in the spin model is 
exactly analogous to the sourceless 't Hooft loop in gauge theories (for a definition see |§). In 
this way the confinement problem can be reduced to Z(2) monopole correlation inequalities. 
Monte Carlo results concerning the relevant monopole correlations will be reported elsewhere 



17| . This gives a unifying picture of confinement in gauge theories and disorder in spin 



models. 

To summarise, we rigorously proved that the presence of a mass gap in the 2d SU(2) 
chiral spin model at low temperature is tantamount to certain vortex correlation inequalities 
in a related SO (3) spin model. These inequalities in turn have been (partly) checked by 
Monte Carlo calculations. It would be worthwile to obtain an analytic proof of the vortex 
correlation inequalities and thus complete the solution of this longstanding problem. Finally, 
it would be also interesting to extend this scheme to other symmetry groups, in particular 
to SU(n). 
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Appendix 



In the Appendix we derive the decomposition of Z e gZ(G c ff — G) into a sum of vortex ex- 
pectations with non-negative coefficients. Using ( Jl8| ) and ([16]) (with the measure du + [U} 
and K(dUi) computed from (p0"[)) this can be written as a double Z(2) integral for two 
independent copies of the uj variables on each plaquette, 



ZZ eS (G eS -G) = fdu + [U] [] fdu p fd^XuM^)^ ex P E^( d *^ + d *^) 



exp K ( d *<^ - d*ui) - exp - ^2 K (d*ui - d*u>i) 
iec 2 iec 2 



We shall make use of the identity 



n ^ - n 9i = s H(fi ± 9i) 



(27) 



(28) 



i=l i=l 

where the summation is on all possible strings of plus and minus signs of length n containing 
an odd number of minus signs. This can be easily verified by induction on n (see ||14j| ). 

Now by applying this identity on the expression in the second line of (p7[) it can be 
written as 



exp K ( d *^ - d *<^) - ex P -K £ ( d *^ ~ d *^) 
iec 2 iec 2 



2L 1 {+-..} iec 2 

We can now expand all the factors in the integrand of fl27D using 



(29) 



f{u) =/(!) !+ co- 



(see equation ([11])) and then make a change of variables 



UJ 



-{uj + uj) 



and 



UJ 



(30) 



(31) 



on each plaquette. This will result in an expression containing factors of the form (uj + ) np {uJ~) n 
on each plaquette and the rest of the integrand will be independent of the a/s. Upon inte- 
grating out the a>'s and cD's each plaquette contributes a factor 



duj n ldu) p (uj+) np {u). 



(32) 
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where 5 npeven =l if n is even, otherwise. This means that each plaquette can have only an 
even power of one type of u (+ or -) otherwise the corresponding term vanishes. 



In the following we list the expansion of different factors in the integrand of (27) in terms 
of the us. Since we want to verify only the positivity of each term in the expansion of (|27 ) 
we shall sometimes omit positive constant factors and in these cases use "~" instead of "=" . 

The plaquette factors can be written as 



Xu v (drip) e hQp = (9 + + 9 u p ) (cosh h + u p sinh h) 



+ 9 u p )(\ + u p h), 



(33) 



where 9^ are projection operators onto states without/with (+/-) a vortex at p, h — tanhh 
and we have omitted the trivial constant factor cosh/i. After the change of variables (|3l|) 
this becomes 



(9 P + to+) u + + (6- - hff})u; + 



-{9; + hB-){u, 
The links not belonging to Ci can be expanded as 



+^2 



+ m-m;)(<o;) 2 (34) 



D Kd*LJi+Kd*u>i 



(1 + Kd k u l )(l + Kd k u l ) 



A-(d*w, + d*tD,) + ^-^(d^,+d*wj) 2 + - — — (d*^ - d*u),) 2 



(35) 



4 v • " 4 

To rewrite this in terms of u^ and u l2 , the uj ± variables on the two plaquettes sharing the 
link I, we can apply the following identies that can be easily checked by using the definition 
of the o^'s; 



LO\U)2 ± OJ1UJ2 = 2(iU^Ujf + UJ 1 UJ2) 



K + ) 2 (^) 2 + fay(u$) 



+ ... 



(36) 
(37) 



The ellipses in the second identity means that terms that integrate to zero because they 
contain both uj + and u~ on the same plaquette have been omitted. Now the final form of 
the expansion for links not contained in C2 is 



2X(^V+ + ^ 2 ) + (l+K 2 )[K) 2 (^) 2 + (cu a ) 2 (^) 



+ (l-K 2 



<4)Vh) 2 + (<*>. 



71YW2? 



(a 



32 



Finally for links contained in C2 it is also straightforward to carry out the same type of 
expansion. These factors can be of two different types depending on the corresponding sign 
in the r.h.s. of ( [29|) carried by the link in question. If the sign is positive we obtain 



e Jf(d*Wi-d*wj) _|_ e -x(d*w;-d*<i>i) _ _ J{ 2 

+2(1 + K 2 



K) 2 (<4) a + (m) 2 (^) 2 



For links with a minus sign the result is 

e X(d^-d^) _ ^(d^-d*^) = + ^f)) 

We are now ready to construct the "diagrams" of our expansion by choosing one term from 
each plaquette and link expansion in all possible ways consistent with the rule that eventually 
each plaquette has to carry an even power of either uo + or uo~ but not both. This constraint 
implies that the plaquettes carrying (9~ ± h9 + )u ± are pairwise connected with stacks of 
links carrying an odd power of the u^'s. These links always come with an additional factor 
~ K while links carrying an even power of the u/s contain 1 ± K 2 . Recall that K = tanhi^ 
is small at low temperature which means that a stack of odd tu-power links costs an energy 
proportional to its length. Moreover the plaquettes that are connected by these stacks carry 
vortices (9~) with a small mixing of the no-vortex state (~ h9 + ). We can thus recover the 
familiar string vortex picture in this representation, the only subtlety being that the effective 
magnetic field gives a small mixing between 9 + and 9~ . 

After integrating out the u variables, up to positive constant factors, each diagram has 
the form 

D = fll(0; ± W+) n_(^ + ± M;)) L , (41) 
peA p^A 

where Qi means expectation with respect to the SO(3) measure di/+[77] and A is a subset of 
plaquettes containing an even number of plaquette 
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